kUAMMMUU 


_  UNCLASSIFIED _ 

security  classification  of  this  rage 


la.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED _ 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 

2b.  OECLASSIFICATI  ON /DOWNGRADING  SCHEDULE 

4  PERFORMING  ORGANIZATION  REPORT  NUMBER1SI 


AD-A162  722 

REPORT  DOCUMENTATION  PAGE 

1  lb.  RESTRICTIVE  MARKINGS 


3.  DISTRIBUTION/AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
unlimited. 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

AFOSR-TR*  r:  120 


6a  NAME  OF  PERFORMING  ORGANIZATION  Sb.  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 

(If  applicable) 

Cfgco  nnivorcit-v  nf  Npu  Ynrlf  Air  Force  Office  of  Scientific  Research 


State  University  of  New  York 

6c.  ADDRESS  (City,  Slate  and  ZIP  Code) 

Applied  Math.  &  Statistics  Dept. 
Stony  Brook,  N:  Y.  11794 


7b.  ADDRESS  (City,  State  and  ZIP  Code) 

Directorate  of  Mathematical  &  Information 
Sciences,  Bolling  AFB  DC  20332 


Sa.  NAME  OF  FUNDING/SPONSORING 
ORGANIZATION 

_ AFOSR _ 

8c.  ADDRESS  (City,  Stale  and  ZIP  Code) 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 

NM _ AFQSR-84-0383 _ 

10.  SOURCE  OF  FUNDING  NOS. 


Bolling  AFB  DC  20332 

11.  TITLE  (Include  Security  CUueification) 

NUMERICAL  METHODS  FOR  DIFFERENTIAL  EQUATIONS 

12.  PERSONAL  AUTHOR(S) 

Dr.  Reginald  P.  Tewarson 


element  no. 


L 

:>  v-v 


A 

.-■-Vi-.] 


L-.v.y. 

A'  ./t 


PROJECT 

TASK 

WORK  UNIT 

NO. 

NO. 

NO. 

2304 

A3 

13a  TYPE  OF  REPORT 


fl3b.  TIME  COVERED 


14.  DATE  OF  REPORT  (Yr  .  Mo  ,  Day)  I  15.  PAGE  COUNT 


16.  SUPPLEMENTARY  NOTATION 


■ll)  COSATI  CODES  [j 

FIELD 

GROUP 

SUB  GR. 

iJ 

| 

mm J 

18  SUBJECT  TERMS  (Continue  on  reverte  if  neceuory  and  identify  by  block  number) 


_ I -  'Computer  code,'  variable  orderf  finite  difference  scheme,^. 

^  ABSTRACT  (Con  unue  on  reverie  if  neccuary  and  identity  by  blue*!  number* 

inuring  thepast  year  the  ®j.nvesLigators  have  been  able  to  develop  a  computer  code  which  has 
[Turned  out  to  be  quite  cofipetitive  with  a  well  established  code.  The  new  approach 
^nplements  a  variable  order  finite  difference  scheme  which  does  not  require  derivatives 
of  the  given  function  and  which  uses  no  informaion  outside  a  subinterval  to  approximate 
the  given  system  in  that  subinterval.  Three  papers  have  been  published  as  a  result  of  this 
effort,  with  the  following  titles,  'An  adaptive  boundary  value  Runge-Kutta  solver  for 
first  order  boundary  value  problems',  'on  the  solution  of  sparse  non-linear  ewuations  and 
some  applications'  adn  1  A  quasi-Newton  method  with  sparse  triple  factorization'. 

Four  additional  papers  are  in  press.  .  ✓  - ",  ^ 


V-f  ,-~t£ 

3  OEC  3  i085 


20.  DISTRIBUTION/AVAILABILITY  OF  ABSTRACT 
UNCLASSIFIED/UNLIMITED  C  SAME  AS  RPT.  □  OTIC  USERS  □ 
22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

JOHN  P,  THOMAS,  Jr.,  Capt,  USAF 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

I  22b  TELEPHONE  NUMBER  22c  OFFICE  SYMBOL 

I  (Include  Area  Code) 

(202)  767-  3026  NM 


DD  FORM  1473,  83  APR 


EDITION  OF  1  JAN  73  IS  OBSOLETE. 


mi 

'-I 


'  V  *. 


l*  *  .  a  ^ 
* 


security  classification  of  this  page 


.  1  . 


nrrn  rr>ilri¥ 


Annual  Report 


Grant  No:  AFOSR-84-0383 


SUMMARY 


The  overall  long  term  goals  of  this  research  project  are: 


to  develop  efficient  computational  algorithms  and  new,  readily 


transportable,  competitive  general  purpose  computer  codes  for 


the  numerical  solution  of  two  point  boundary  value  differential 


equations;  to  use  these  techniques  to  develop  adaptive  codes 


stiff  and  non-stiff  initial  value  problems  and  the  method 


of  j.  .  for  the  solution  of  partial  differential  equations; 


to  test  he  codes  on  problems  arising  in  diverse  application 


areas,  e.g.,  Physics,  Chemistry,  Fluid  Mechanics,  Combustion 


and  Biomatheraati. 


The  short  term  cal  is  to  develop  and  implement  variable 


order  finite  difference  schemes  which  do  not  require 


derivatives  of  the  given  function  and  use  no  information 


itside  a  subinterval  to  approximate  the  given  system  of 


difi.  -cential  equations  in  that  sub.,  terval  and  use  a  variation 


of  defe  '•ed  correction  for  the  soluti<_  of  the  corresponding 


non-linear  .  1qebraic  equations. 


During  the  oast  year  the  investigators  hare  been  able  to 


develop  a  computer  c.-’de  which,  on  the  basis  of  some  preliminary 


experiments,  has  turned  out  to  be  quite  competitive  with  a  well 


established  code.  Extensive  test  are  planned  during  the  coming 


year.  The  algorithmic  development  phase  of  the  research  has  led 
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to  several  useful  results.  An  eighth  order  method,  which  has 
second,  fourth  and  sixth  order  methods  embedded  in  it  has  been 
developed.  The  method  possesses  some  of  the  best  features  of 
implicit  Runge-Kutta  and  Gap  schemes.  A  multiderivative 
generalization  of  the  above  schemes  has  also  been  realized.  A 
sparse  factorization  for  quasi-Newton  type  methods  has  been 
obtained.  Algorithms  for  solving  combined  systems  of  linear  and 
non-linear  algebraic  equations  and  matrix  splitting  have  also 
been  developed. 

Some  of  our  results  have  been  used  for  the  computer 
simulation  of  kidney  function  at  the  Cornell  Medical  School  in 
a  collaborative  research  in  which  the  Principal  Investigator  is 
involved . 

REPORT 

Finite  difference  methods  which  combine  features  of  both 
Runge-Kutta  processes  and  Gap  schemes  have  been  developed  [1]. 
These  methods  are  suitable  for  use  in  adaptive  codes  for  the 
solution  of  first  order  differential  equations  with  two-point 
boundary  conditions.  Success  has  been  achieved  in  finding  the 
order  conditions  for  these  methods  and  also  some  techniques 
that  can  be  used  to  reduce  the  order  number.  An  eighth  order 
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method  which  has  second,  fourth  and  sixth  order  methods 
embedded  in  it  has  been  developed.  This  eighth  order  method  as 
well  as  embedded  methods  are  A-stable.  A  technique  has  been 
proposed  which  exploits  the  natural  embedding  of  methods  for 
the  estimation  of  the  error  in  the  numerical  solution.  This 
natural  embedding  of  the  methods  is  also  utilized  to  estimate 
the  local  truncation  error.  In  turn,  these  truncation  error 
estimates  are  used  to  generate  asymptotically  equidistributing 
meshes  using  adaptive  mesh  placement  procedures.  A  quasi- 
Newton  method  is  described  in  this  paper  for  solving 
discretized  systems  arising  from  the  use  of  proposed  methods. 
It  has  been  possible  to  developed  a  preliminary  general  purpose 
adaptive  code  based  on  these  methods. 

A  general  description  of  sparsity  based  algorithms  that 
use  structural  information  in  handling  the  model  equations  is 
given  in  [2].  The  primary  focus  is  on  splitting  of  equations 
and  variables  into  two  subsets:  one  large  but  relatively  easy 
to  solve,  called  the  non-basic  subset,  and  the  other  small  but 
difficult  to  handle,  called  the  basic  subset.  The  non-basic 
subset  is  then  solved  more  often  than  the  difficult  basic  set. 
The  computer  storage  is  that  required  by  the  basic  set. 
Applications  in  flow  networks  and  energy-economy  models  are 
briefly  described. 
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A  new  quasi-Newton  method  for  sparse  matrices  arising  in 
various  application  areas  is  presented  in  [3].  In  this  work 
sparse  triple  factorization  is  used  to  develop  a  rank  two 
update . 

In  [4],  multi-derivative  generalizations  of  the  processes 
given  in  [1]  are  considered.  These  use  fewer  derivatives  to 
yield  results  analogous  to  those  for  Gap  schemes.  Obviously, 
these  processes  instead  of  Gap  schemes  should  be  utilized  if 
the  higher  order  analytic  derivatives  are  not  readily 
available . 

An  efficient  method  for  handling  systems  with  linear  and 
non-linear  subsystems  is  given  in  [5].  In  this  method 
iterations  are  needed  only  for  the  solution  of  non-linear 
subsystems . 

Cubic  and  quintic  splines  have  been  utilized  in  some  of 
the  methods  for  solving  differential  equations  developed  by  the 
investigators  and  others  in  the  past  .  One  of  the  crucial 
steps  in  these  techniques  involves  the  solution  of  band 
systems.  A  very  efficient  general  purpose  scheme  for  handling 
such  band  systems  [6]  has  been  developed.  This  scheme  is  also 
very  useful  for  block  banded  systems  arising  in  a  variety  of 
applications . 
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A  method  of  order  eight  that  uses  cubic  splines  on  quintic 
splines  for  solving  differential  equations  is  described  in 
[7], 

A  class  of  processes  [8]  have  been  found  for  solving 
second  order  boundary  value  problems  of  the  form  y''  -  f(x,y), 
y(a)«A,  y(b)-B.  These  do  not  require  derivatives  of  f,  and 
result  in  nonlinear  equations  with  blocked  tridiagonal 
Jacobians.  Specialized  techniques  have  been  developed  to 
derive  the  necessary  formulas.  The  computation  of  the  Jacobians 
for  the  discretized  systems  arising  from  the  use  of  these 
formulas  is  expensive.  This  makes  the  ordinary  Newton 
iterations  impractical.  This  paper  contains  various  techniques 
for  approximating  and  splitting  the  Jacobians.  These 
techniques  are  computationally  inexpensive  and  lead  to 
excellent  convergence. 

The  Principal  Investigator  has  given  invited  papers  based 
on  some  of  the  research  resulting  from  this  grant  at  two 
international  conferences  [9,10].  One  of  the  problems  in  the 
solution  of  large  systems  of  non-linear  algebraic  equations 
arising  in  the  numerical  solution  of  differential  equations  is 
to  be  able  to  choose  the  starting  approximate  solution  in  the 
domain  of  attraction  of  the  desired  solution.  Continuation 
type  methods  are  used  to  enlarge  this  domain  of  attraction.  An 
algorithm  (based  on  Moore-Penrose  generalized  inverse)  to 
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handle  the  turning  points  in  continuation  methods  has  been 
developed. 
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